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Abstract 

From the point of view of stochastic analysis the Caputo and Riemann-Liouville 
derivatives of order a € (0, 2) can be viewed as (regularized) generators of stable 
Levy motions interrupted on crossing a boundary. This interpretation naturally 
suggests fully mixed, two-sided or even multidimensional generalizations of these 
derivatives, as well as a probabilistic approach to the analysis of the related equa¬ 
tions. These extensions are introduced and some well-posedness results are obtained 
that generalize, simplify and unify lots of known facts. This probabilistic analysis 
leads one to study a class of Markov processes that can be constructed from any 
given Markov process in R'^ by blocking (or interrupting) the jumps that attempt 
to cross certain closed set of ’check-points’. 
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1 Introduction 

For a general background in fractional calculus and fractional equations we refer to books 
0, ra. P!. H. see also survey [M], for the crucial link with CTRW to |32], [16] and 
m and for the numerous applications in natural science to [33] and [M] . 

The aim of this paper is to present a systematic treatment of a class of equations 
that include fractional derivatives as very particular cases. Unlike the mostly analytic 
studies of fractional differential equations (see the reference above), the present treatment 
and the corresponding far reaching extensions of fractional derivatives are based on a 
probabilistic point of view. This link with probability provides a powerful tool for the 
study of fractional equations. It is also worth mentioning the recent activity on proving 
probabilistic interpretation of solutions by analytic methods, see e.g. [TT], while our 
approach provides such an interpretation as a starting point. 
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From this point of view the basic Caputo and Riemann-Liouville (RL) derivatives 
of order a G (0, 2) can be viewed as (regularized) generators of stable Levy motions 
interrupted on crossing a boundary. This interpretation naturally suggests fully mixed, 
two-sided or even multidimensional generalizations of these derivatives, as well as a proba¬ 
bilistic approach to the analysis of the related equations. These extensions are introduced 
leading to well-posedness results that generalize, simplify and unify lots of known facts. 
Some explicit solutions are also obtained. The corresponding probabilistic analysis leads 
one to study an interesting general class of Markov processes that can be constructed from 
any given Markov process in by blocking (or interrupting) the jumps that attempt 
to cross certain closed set of check-points. This analysis is only initiated in the present 
work. Further development, as well as the application of this technique to the study of 
fractional in time and space diffusion equations and to the convergence of CTRW, will be 
discussed in separate publications. 

The paper is organized as follows. In the next preliminary section we introduce in a 
convenient form the main objects of fractional calculus, the Caputo and RL derivatives. 

In Section [3] we explain in detail the probabilistic meaning of these derivatives and 
their natural place in stochastic analysis leading, on the one hand side, to far reaching 
generalizations, and on the other hand, to a unihed treatment of various equations by 
powerful tools of stochastic analysis. We distinguish the cases of /5 G (0,1) and [3 G (1, 2), 
because in the hrst case the Caputo derivative has a direct probabilistic interpretation 
and in the second an additional regularization is needed. We also treat separately a multi¬ 
dimensional extension, as it includes one additional ingredient, the projection of a jump 
on a boundary. 

The hnal Section H] initiates the rigorous theory of the equations and processes intro¬ 
duced above by providing some basic examples (not at all exhaustive) of well-posedness 
results and explicit solutions that can be obtained by these tools, with main attention 
restricted to the analogs of the derivatives of order [3 <1. 

Appendix derives some equivalent versions of basic fractional derivatives mentioned 
in the next section without proof. These calculations should be obvious for specialists in 
fractional calculus and are given here for completeness. 

We shall denote by 1 m the indicator function of a set M. 


2 Preliminaries: classical fractional derivatives 

For convenience we recall here the basic dehnitions of fractional derivatives and their 
equivalent representations htting our purposes. 

We shall repeatedly use the Euler identity F(a:) = {x — l)F(x — 1) for the Euler 
Gamma-function whenever appropriate without mentioning it. 

Due to the formula for the iterated Riemann integral 

w = feVT)! ~ (1) 

it is natural to extend it analytically, if a; > a, to complex n with positive real part, leading 
to the following dehnition of the (right) fractional or Riemann-Liouville (RL) integral of 
order f3 (with positive real part): 

tf~^f{t)dt. (2) 
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Notice that for x < a, we have negative numbers {x — t) in power n — 1, so that the 
corresponding extension to complex (or even real) n leading to the so-called left RL 
integrals have some subtleties to be discussed later. 

Noting that the derivation is the inverse operation to usual integration, the above 
dehnition suggests two notions of fractional derivative, the so-called RL (right) derivatives 
of order > 0, /3 ^ N; 


dU{^) = 


dx^- 


_ jn-p 
) o.-\- 


fix) = 


d” 


r(n — 13) dx'"- 


{x — tY ^ ^f{t)dt, X > a, (3) 


and the so-called Caputo (right) derivative of order /3 > 0, /3 ^ N: 


dLj(x) = 4"-'’ 


dx' 


:f 


[X = 


r(n-/3) 


[X 


i) 


n-P-l 


—/ 
dt^ 


{t)dt, X > a (4) 


where n is the maximal integer that is strictly less than (3 + 1. 

Straightforward integration by parts (see Appendix) show that, for smooth enough / 
and /9 G (0,1), a; > a. 




1 r- fix - z) - nx) f(x) 

T{-l3) Jq z^+I^ r(l-/3)(a;-a)^’ 


(5) 


implying 


DtJix) 


1 r ° fix -z) - f{x) f{x) - /(g) 

r(-/3)io ^ T{l-/3){x-aY^ 


D^a+Jix) = - fia)]{x) = Df+/(a;) 


/(q) 

T{1- /3)\x-a\f^' 


( 6 ) 

(7) 


In particular it follows that for smooth bounded integrable functions, the right RL 
and Caputo derivatives coincide for a = —oo, (3 G (0,1), and one dehnes the fractional 
derivative in generator form as their common value: 




f{x) = Dt^^f{x) = Df^^J{x) 


1 fjx-z) - fix) 

r(-/9) Jo 


( 8 ) 


Analogously (see Appendix for detail), for smooth enough / and (3 G (1,2), x > a, 
one hnds that 


L>f+/(x) = 


r(-/9) 


f{x -z) - f{x) -1- f{x)z , fix){x - a) ^ , I3f'ix){x - af ^ 




D^a+Jix) = 


r(-/9) 


r(i-/3) 

f{x -z) - f{x) + f{x)z 
R+3 


r(2-/3) 


dz 


(fix) - f(a))(x - a)_p_ (fsfix ) - f'{a))(x - a)‘ f 


r(l-/?) 


(9) 


( 10 ) 


SO that 


L>f+*/(x) = L>f+[/-/(a)-/'(a)(.-a)](x) = D^^f{x) 


T{2-/3) 

/(a)(x-a)"^ /'(a)(x - a)^"^ 


r(i-/9) 


T{2-/3) 


( 11 ) 
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Again for smooth bounded integrable functions, the right RL and Caputo derivatives 
coincide for a = —oo, (3 G (1,2), and one dehnes the fractional derivative in generator 
form as their common value: 


dx^ 


f{x) = = C-oo+i/la:) 


1 /’” f(x- z) - f(x) + f'(x)z 

n-p) I 


( 12 ) 


Turning to the left derivative notice that for x < a formula ([I]) rewrites as 


4 ”/(^) = l‘(t - xr-'f{t)dt, ( 13 ) 

suggesting several possible normalizations for the analytic continuation in n and the cor¬ 
responding inversions (fractional derivatives). We are interested only in the derivatives 
of order less than 2. For a unihed probabilistic interpretation of these derivatives it is 
convenient to choose the left versions of ([3]), (jl]) as follows. For (3 G (0,1): 


D^_f{x) = j ^ 

Da-J{x) = dt, X <a] (15) 

for/? G (1,2): 

D^_f{x) = 13 -^ dx^ I dt, x<a, (16) 

D^Jix) = y{^2- 13) J dt, X <a. (17) 

When /3 G (0,1) and x < a, similar calculations as for the right derivative (see fIlOQp l 
lead to the following analogs of (j5]), ([6]): 


Dtnx) 


1 r-* f{x + z)- fix) f(x) 

r(-/?)io T{1 - fS){a - xY' 


(18) 


DLtf(x) 

implying 


1 f{x + z)-f{x) fix) - /(g) 

Ti-^)Jo F(l-/3)(a-a;)/5’ 


F>f_*/(x) = F>f_[/ - /(a)](x) = D^fix) 


/(a) 

T{1- I3)ia-x)f^' 


(19) 


( 20 ) 


When /3 G (1, 2), X < a, one obtains 


DUix) = 


r(-)9) 


fix + z)- fix) - ^ fix)ia - x) ^ (3f'ix)ia-x)^ ^ 


Dtjix) = 


z^+P F(1 - /3) 

r~^ fix + z)- fix) - f'ix)z 


Ti2-(3) 


( 21 ) 


r(-/3) 


'0 




1+/3 


dz 
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( 22 ) 


{fix) - f{a)){a-x) ^ {f3f{x)-f{a)){a-x) 


+ 






so that 


D^_J{x) = D^_[f-f{a)-f{a){.-a)]{x) = D^_f{x) 


— 


r(2-/3) 


f{a){a — x) ^ f'{a){a — x)^ ^ 


r(i-/9) 


r(2-/3) 


(23) 

For smooth bounded integrable functions the left fractional derivatives in generator 
form become 


^/3 1 

f(x) = = Di_J(x) = 


d{—x)^ 


r(-/9) 


(24) 


_JP_ 

d{—x)f^ 


fix) = D^-fix) 


DL-Jix) 


1 /‘°° fix + z)- f{x) 

r(-/d) Jo 


f'ix)z 


dz, (25) 


for fd G (0,1) and (d G (1, 2) respectively. 

It is straightforward to see that the pairs of operators (IH]), ll2^ and flT2l) . fl25ll are dual 
in the sense that 




d^ \ 
d{—x)f^^) 


for fd G (0,1)U(1, 2) and sufficiently regular functions /, g, where the pairing (/, g) denotes 
of course the usual L^-product: if,g) = f fix)g(x)dx. This fact also justifies the notation 
d^/d{—x)^, since for = 1 the operators d/dx and —d/dx = d/d{—x) are dual. 


3 Probabilistic interpretation and extensions of RL 
and Caputo derivatives 


3.1 The case ^ £ (0, 1), (i = 1 

It is well known that the operators —d^/d{—x)^ and d^/d{—x)^ from fl2T)) and (|2^ respec¬ 
tively are the generators of stable Levy motions without negative jumps, the (annoying) 
discrepancy in the sign reflects the fact that r(— /3) < 0 for G (0,1) and r(— /?) > 0 for 
fd G (1, 2). In particular, for fd G (0,1), the corresponding Levy process is a stable subordi- 
nator (an increasing process). Similarly, the operators —d^/dx^ and for fd G (0,1) 

and fd G (1, 2) respectively generate stable Levy motions without positive jumps, which 
can be obtained by inversion of the processes generate by —d^/d(—x)^ and d^/d(—x)^. 

Let us now look specifically at the decreasing process X{t) generated hj A = —d^/dx^, 
fd G (0,1). Let us modify it by forbidding it (interrupting on an attempt) to cross a 
boundary x = a with an a G R, that is, all jumps aimed to land to the left of the 
chosen barrier-point a are forced to land exactly at a. Analytically, this procedure means 
changing the generator A to 


Aa+*fix) 


1 r-^ fix-z)-fix) 1 r fja) - fix) 

r(-/d) Jo z^+^ " r(-/?) 


X > a, 


which rewrites as 


^a+Jix) 


1 r “ fix -z) - fix) fja) - /(x) 

ri-fd)Jo z^+f^ T{l-/d)ix-a)0 


X > a. 


(26) 
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In this expression we recognize the Caputo derivative (jl]), with inverted sign. Killing the 
process at the boundary-point x = a means analytically to set /(a) = 0, in which case 
fl26|) turns to RL fractional derivative (with inverted sign). 

We conclude that the transition from the free derivative /dx^ to the Caputo right 
derivative at a is a particular case of the procedure of interrupting a decreasing process 
on crossing a boundary. Namely, let an operator 


Af{x)= {f{x-y)-f{x))u{x,dy) 


with a kernel ^{x ,.) on (0, oo) such that 


sup / min(l, \y\)u{x, dy) < oo 

X Jo 


(27) 


generate a decreasing Feller process. Then the corresponding process interrupted (and 
hence stopped) at a boundary x = a has the generator 

nx—a nOQ 

Aa+^f{x)= {f{x-y)-f{x))iy{x,dy) + {f{a)-f{x)) i^{x,dy), x > a. (28) 

Jo J x—a 

Similarly, the transition from the free derivative d^/d{—x)^ to the Caputo left deriva¬ 
tive at a is a particular case of the procedure of interrupting an increasing process on 
crossing a boundary, that is, the transition from a process on R generated by 

poo 

Af{x)= {f{x + y) - f{x))iy{x,dy) 

Jo 

to the process on the interval (—cx), a] generated by 

pa—x poo 

Aa-i.f{x)= {f{x+ y) - f{x))iy{x,dy) + {f{a) - f{x)) i^{x,dy), x < a. (29) 

Jo Ja—x 

From this point of view, the natural extension of this procedure to the general processes 
of bounded variation generated by the operators 


Af{x) ='y{x)f'{x) + / {f{x + y) - f{x))iy{x,dy) 

J —OO 

with a kernel v{x ,.) on R \ {0} such that 


(30) 


sup / min(l, \y\)v{x^dy) < oo 

X J-oo 


(31) 


is a transition to the process in a given interval [a, h] with arbitrary —cx) < a < b < oo 
with jumps interrupted on an attempt to cross the boundary (i.e. all jumps aiming to 
land outside [a,b] are forced to land on its nearest point), that is the process generated 
by the operator 


f*b—x 


A[a,b]i,f{x) = 'y{x)f{x) + / {f{x + y)- f{x))u{x, dy) 
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poo pa—x 

+ {f{b)-fix)) i/{x,dy) + {f{a) - f{x)) v{x,dy), x e {a,b). (32) 

J b—x J —oo 

To deal with fractional derivatives we are mostly interested in the motion inside [a, b]. 
However, for the analysis it is convenient to be able to have the corresponding process 
extended to all R. This can be done in various ways, though two natural approaches can 
be distinguished. In the hrst one all jumps are supposed to be restricted to jump in the 
direction of D only, for instance one can choose the extension of to all x is given 

by the expression 


A[a,b]*f{x) = 'y{x)f{x) + la;<b / {f[{x + y) Ab]- f{x))u{x, dy) 

Jo 

+ lx>a f {f[{x + y)y a]-f{x))u{x,dy). (33) 

J —OO 

In the second approach we stick to the idea of interruption on crossing the boundary, so 
that D and its complement are treated symmetrically. To present this approach in a proper 
generality assume a hnite set R = {6i < ■ • • < bk) or a countable set B = {bi,i G N}, 
with bi < for any i, is chosen. For instance, in the above setting R is a two-point 
set R = {a, b}. For any a; G R let us now dehne b^{x) and 6_(x) as the nearest point of 
R to the right and to the left of x respectively {x excluded in both cases even if a; G R). 
Then the modihcation of the process on R generated by fl5U]l . with jumps interrupted on 
crossing R (think of R as a set of road blocks or check points placed to control the free 
motion given by A) can be specihed by the generator 


^B*f{x)='y{x)f{x)+ {f{[{x + y) Ab+{x)]\/b_{x))-f{x))iy{x,dy). (34) 

J —OO 

Clearly for x G (6i,6i+i) this generator coincides with fl32l) . where a = bi,b = fej+i. In 
order to study the fractional differential equations in intervals we can choose to work with 
the processes given by (|33ll or (l3T)l . as their behavior until they reach the boundary dD 
is identical. For more general domains D these different extensions can lead of course to 
different problems. 

In [18], [19] the author referred to generators of type (130D as to the generators of 
order at most one aiming to stress that they can be considered as fully mixed fractional 
derivatives of order not exceeding one. 

If we kill the process generated by A at the boundary, the generator (15^ turns to 



{f{x + y)- f{x))u{x, dy) - f{x) 


v{x,dy) 


u{x,dy) 


_J b—x 


, (35) 


for X G (a, 6), which represents the corresponding extension of RL fractional derivative. 

Therefore, from probabilistic point of view, the natural extension of the basic linear 
fractional equation 

Da+J{x) = -A/(a;), x > a, (36) 

with initial condition fa = /(a) and A > 0 given, is the problem of hnding / on [a, oo) 
satisfying 

Aa+J{x) = Xf{x), f{a) = fa, 
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with Aa+-k given by fl28|) . that is the equation 



y) - dy) + (/(a) - f{x)) 


u{x,dy) = Xf{x), 


/(a) = U (37) 


This equation includes the extensions of (I5B]) with various mixed fractional derivatives, 
that is, the problem 


'^UjD^i^f{x) = -Xf{x), f{a) = fa, (38) 

j 

with some hnite collection of numbers Uj > 0 and fdj G (0,1) or even more exotic versions 
with Uj or fdj being functions of x. Mixed derivatives and related fractional equations are 
actively studied recently by analytical methods, see e. g. [7], [22], [30], [I2] . 

Moreover, it is now natural to formulate the two-sided version of fl36|) as the equation 

+ Dljix) = f{a) = /„, f(b) = f, (39) 

on the interval x G [a, b], which represents the simplest case of a more general equation 

Aia,b]*f{x) = Xf{x), f{a) = fa, f{b) = fb (40) 

with given by (l33D . In particular, the mixed-derivatives extension of fl39H is the 

problem 


k k 

- '^ujjD^aUfi^) -'^IjDbfJix) = Xf{x), f{a) = fa, f{b) = fb, (41) 
i=i i=i 

with some collection of numbers (or, more generally, functions) Uj > 0, jj > 0, fdj G (0,1) 
and A > 0 extending for instance the problem considered in HD]. The well-posedness for 
this problem is covered by Theorem 14.21 below. 

As solutions to the equation /!„+*/(x) = g{x) with a given g and initial condition 
f{a) = a can be given in terms of fractional integrals, the two-sided analog of fractional 
integral becomes the solution to the problem 

Da+*f{x) + D^_J{x) = -g{x), f{a) = fa, f{b) = fb, (42) 

with given fa, fb and g on [a, b], which again represents the simplest case of a more general 
problem 

X^[a,b]Jix) = g{x), f{a) = fa, f{b) = fb, (43) 

or with the linear term included 

A[a,b]*f{x) = Xf{x) + g{x), /(a) = fa, fib) = fb, (44) 

suggesting a new class of extensions of the notion of fractional integral, which is alternative 
to a more classical one, see [15], which is based on the variations of special functions used 
as the kernels (which is natural from an analytic perspective). 




Remark 1. Notice that equation can be rewritten in terms of RL derivative. Namely, 
let v{x) be a smooth function on [a,b] with v{a) = fa,v{b) = fb. Then for the function 
(p = f — V problem rewrites as 

-4[a,b]*0 = \(t) + g- A\a,b\i.v + \v, 0(a) = 0(6) = 0, (45) 

which is equivalent to 

A[a,b](p = X(p + g, 0(a) = 0(6) = 0, g = g + {X- A[a,b]*)v. (46) 

Particular examples of boundary value problems with fractional derivatives are actively 
studied, mostly by analytical techniques, see e. g. [13], [H] or [26] and references therein, 
for distributed or mixed derivatives see also ra. ni and [35] . 


3.2 The case {3 £ (0,1), d > 1 

Let us now turn to the multidimensional extension of this interruption procedure. The 
analog of RL derivative arising from a process in and a domain D C R*^ is the 
generator of the process killed on leaving D. For Caputo version this is more subtle, as 
we have to specify a point where a jump crosses the boundary. Below we choose the most 
natural model assuming that a trajectory of a jump follows shortest path (a straight line 
in Euclidean case). Suppose ^4 is a generator of a Feller process Xt{x) in R'^ with the 
generator of type 

Af{x) = {^{x),V)f{x) + [ {f{x + y) - f{x))u{x,dy) (47) 

with a kernel u{x,.) on R'^ \ {0} such that 


sup / min(l, ||/|)z/(a;, d?/) < oo, (48) 

a; JRd 

that is, in the terminology of [18], [19], a generator or order at most one. 

Let D be an open convex subset of R'^ with boundary dD and closure D. For x G R'^ 
and a unit vector e let = {x + Ae, A > 0} be the ray drawn from x in the direction e. 
For X G R'^, let 

D{x) = {y eR'^ : Lx,y/|y| n i) 7^ 0}. 

In particular, D{x) = D for all x G iA. Furthermore, for y G D{x), let 

A(x, y/\y\) = max{i? > 0 : x + Ry/\y\ G D}. 

Let us introduce the restriction function R£,{x,y), x G R*^, y G D{x), by the formula 


x + y, if ||/| < A(x, 2 //||/|) 

X + X{x,y/\y\)y/\y\, ii\y\ > X{x,y/\y\) 


(49) 


The process Xt{x) with jumps interrupted on crossing dD when jumping from inside D 
and forced to jump in the direction of D when jumping from outside of D can be defined 
by the generator 


AD*f{x) = ( 7 (x), V)/(x) + f [fiRnix, y)) - /(x)]z/(x, dy), (50) 

Jd{x) 
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which represents a multidimensional extension of the Caputo boundary operator in D 
arising from fl33|) . 

To define a multidimensional analog of fl3T|) let us assume more generally that D is 
an arbitrary open subset of R'^, so that B = \ i) is closed. Let us define 

X{x, y/\y\) = min{R > 0 : a: + Ry/\y\ G B} 

for any y, with the convention that X{x,y/\y\) = cx) if the ray does not intersect 

B at all, and let the projection-on-the-boundary function RB{x,y) be defined for all 
x,y E R'^ by the same formula (09]), but with A instead of A. Then the analog of fIMll . 
that is the modification of the process on R'^ generated by (071) . obtained by interrupting 
jumps on an attempt to cross R, is specified by the generator 

= ( 7 ( 3 ^), V)/(a:) + /" [f{RD{x,y)) - f{x)]u{x,dy). (51) 

aRd 

For a convex domain D the operators fl50|) and (l5T|) with B = dD and D = R'^ \ B 
coincide for x E D, so that when one is interested in the random motion inside D one can 
work with either of the processes generated by fl50|l or (l5T|) and stopped on the boundary 
of R. 

As pointed out above, the process Xt{x) killed on the boundary, which represents 
a multidimensional extension of the RL boundary operator in D arising from (06|) . is 
specified by the generator 


^Df{x)=[ [f{x+ y)l^+y^D - f{x)]jy{x,dy), x E D. (52) 

jRd 

Remark 2. In dimension d = 1 the projection z — )■ R(^a,b){x,z) = [z A a] \/ b of the 
real numbers to the interval [a, 6 ] (used in fl33|) ) does not depend on x and is clearly the 
most natural one. In higher dimensions one can imagine several reasonable extensions. 
Our choice used above was meant to preserve the direction of a jump. Another reasonable 
choice could be the definition of the projection Rd{x, z) as the point on D nearest to z (both 
choices coincide in d = 1). This would lead to a different multi-dimensional extension of 
the Caputo derivative. 

Assuming for simplicity that the kernel u has a density, iy(x;y), with respect to 
Lebesgue measure, consider the following three basic examples. 

If R is a half space 

D = Db = {{xi,X 2 ) E R'^^^ : Xi <b,X2 E R'^}, (53) 


then 


h — Xf 

Rd{x, y) = Rd{x, y)= (b,X2^ - 1/2 1 , xi < b < xi + yi, 


yi 




Aoi^fix) = / dy2 

JR'i 


dyiu{x;y)[f{x + y) - f{x)] 


b — Xi 


+ / dy 2 dyiiy{x-,y)[f [b,X 2 -{ - 2/2 -/(x)]. 


J b—xi 


yi 


(54) 
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If D is a band 


D = D(a,b) = {{xi,X 2 ) G : a < Xi < 6 ,X 2 e R'^}, 
then, for Xi G (a, b), 

p pb—xi 

ADi,f{x)= dy 2 dyiu{x-,y)[f{x + y) - f{x)] 

J J a—x\ 

+ [ dy 2 [ dyiiy{x;y)[f (b,X 2 + - — —y 2 ^-f{x)] 

jR-i Jb-xi V Dl J 


f*a—xi 


+ 1 dy2 I dyiu{x]y)[f (a,X2 + - — —y2]-f{x)]. 

•/r^ j—oo V 2/l / 

If D is the unit ball D = {x \ \x\ < 1} in R'^, then, for x G -D, 


(55) 


(56) 


Rd{x, y) = Rd{x, y) = x+\{x, y)y, A(x, y) = 

\y\ 


^1,1 I 12 / y 

"’wJ h'M 


AdJ{x) = dy z/(x; y) (lA(a;,.)>i [fix + y) - fix)] + 1a(x,.)<i [fix + A(x, y)y) - f{x)]) . 

(67) 


3.3 The case /? € (1,2) 

Let us now look at the derivatives of order /3 G (1,2). Interrupting the Levy process with 
only negative jumps generated by flT^ on crossing the boundary {x = a} means changing 
its generator to the operator 


Dlj(x) = 


which rewrites as 


fix - z)- fix) + f'ix)z ^^ ^ 1 “ 


/(a) - f{x) + f'{x)z 


yl+/3 


r(-/9) 




dz, 


(58) 


Da+fix) = 


r(-/9) 


f{x - z)- f{x) + f{x)z 


dz 


^ ifjx) - fia))ix-a) y ^ ^fix){x - a)^ ^ 


r(i-)d) r(2-/i) • 

This expression differs from the Caputo derivative ffTOj) by the term containing f'ia): 

f'ia){x - a)^-y 


(59) 


D^^^Jix) = D^^^fix) 


r(2-/l) 


(60) 


which one could expect as operator (na does not have a structure that allows one to 
interpret it as a generator for a Markov process precisely because of this term containing 
f'ia) (recall that r(— /3) > 0 and r(l — /S) < 0 here). 
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In order to see the meaning of this correcting term let us observe that if / G C^[a,oo), 
then, up to terms tending to zero, D^+fix) behaves like 


f\x){x-aY ^ 


+ 


/S 


Lr(i-/3) r(2-/3) 


f\x){x — a)^ ^ 

r(2-/3) ’ 


as a: —>■ a, and thus tends to ±oo if f\a) is positive or negative respectively. Thus 
subtraction of the term containing f'{a) in fl60ll is a regularization of that makes it 
hnite on smooth functions. 

On the other hand, killing the process generated by flHHll at the boundary point x = a 
means setting /(a) = 0 and then fIbOll turns exactly into the Riemann-Liouville derivative 


D^ix) = 


r(-/9) 


f{x -z)- f{x) + ^ f{x){x-a) ^ ^ /3f'{x){x - af ^ 


yl+/3 


r(i-/3) 


r(2-/3) 


(61) 

precisely as in the case (3 G (0,1). 

Extension of this procedure is thus clear. Namely, starting with a Feller process with 
negative jumps, say with the generator 


^f{x)= / [f{x-z)-f{x) +f{x)z]u{x,dz) 


such that 


sup / z‘^u{x, dz) < oo, 

X Jo 


zui^x, dz) = oo 


(62) 


(63) 


one can form the corresponding process with jumps interrupted at a as the process gen¬ 
erated by 

Aa+f{x) = [ [f{x - z) - f{x) + f\x)z]u{x, dz) 

Jo 

POO PCX) 

+ (/(a) -/(a;)) / u{x,dz) + f{x) zu{x,dz). (64) 

Jx—a Jx—a 

Its main term of asymptotics as x —>■ a is 

POO 

f{x) / {z-{x- a))u{x,dz), 

J X—a 

which is unbounded unless f'{a) = 0. Thus the analog of the Caputo fractional derivative 
is obtained by subtracting this ’inhnity’: 


Aa+J{x) = Aa+f{x) - f{a) / (z- (x- a))z/(a, dz). 

J X—a 

Moreover, it implies the following limiting behavior for / G (^^[a, oo): 

lim Aa+*f{x) = 0, 

X—>-a+ 

lim Aa+f{x) = 0, if/'(a) = 0. 

X— 


(65) 


( 66 ) 
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Remark 3. Notice that generator (IM|) has variable coefficients even if underlying process 
was a Levy process. Hence it is not directly clear whether (IMD generates a well defined 
process. This issue will be addressed in the next sections. 

Similarly one defines for Feller processes with positive jumps. More generally, 
for any Feller process on R generated by the operator 

Af{x) = G{x)f{x)+-f{x)f{x) + [ [f{x +z) - f{x) - f{x)zxiz)]iy{x,dz), (67) 

Jn 

where x is some molliher (that traditionally is taken as either l|z|<i or as 1/(1 + z"^)), one 
dehnes the corresponding process with jumps interrupted on crossing an interval [a, b] as 
the process generated by the operator 

+ l{x)f{x) + f [f{x + z)- f{x) - f\x)zx{z)]u{x, dz) 

J a—x 


+ {f{b)-f{x)) / u{x,dz) + {f{a)-f{x)) / u{x,dz)-f{x) / 

Jb—x J—oo J'R\{a—x,b—x) 

for X G (a, b) (when such process is well dehned), or for arbitrary x 


zx{z)u{x,dz) 

( 68 ) 


Aia,b]f{x) = G{x)f'\x)+x{x)f\x)+l^<b / {f[{x + z) Ab]- f{x) - f{x)zx{z)) u{x, dz) 


^x>a 


(/([(a: + 2 ) V o] - f[x) - f'{x)zx(z))v(x, dz). 


(69) 


The corresponding analog of the Caputo derivative is obtained by subtracting the 
singularity at the boundary points, that is, if 


J \z\i'{x, dz) = 00 , I zi'{x, dz) = 00 , 

then, for x G (a, b), 

poo 

AxihfiO = / Hl>-z:)-zx{z)]iz(x,dz)-f'{a) 


(70) 


' b—x 


[(a-x)-zx(z)]iz(x, dz). 

(71) 


The analog of the Riemann-Liouville derivative is obtained from A[a,b]f{x) by setting the 
boundary values of / to zero yielding the operator of the processes generated by A[a,b]f{x), 
but killed at the boundary: 


f*b—x 


7l[a,b]/((E) = G{x)f{x) + x{x)f\x) + / [f{x + z)- f{x) - f\x)zx{z)]u{x, dz) 


f{x) iy{xffiz)-f{x) zx{z)iy{xffiz). 

J'R\{a—x,b—x) J'R\{a—x,b—x) 


(72) 


But what is the relation between equations involving A\^a,b]-k and A[a,bfi The point is 
that if A[a,b]f = 9 on [a,b] for some bounded g, then necessarily f'{a) = f'{h) = 0 (as 
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otherwise Aya,b] would be unbounded for a; —)■ a or x —)■ 6 ) and hence A[a,b]f = A[a,b]*, so 
that at least classical solutions for the problems 

= ^/(^) + 9 {x), f{a) = fa, f{b) = fb (73) 

also solve the problem 

A[a,b]i,f{x) = Xf{x) + g{x), f{a) = fa, f{b) = fb- (74) 

The problem fl73|l can be naturally settled probabilistically. 

Remark 4. The subtraction of singularity in definition CD makes it dependent on repre¬ 
sentation dnn). This is however not very essential for solving the corresponding fractional 
differential eguations, because operator A is not representation dependent. 

Similarly in multidimensional case the analog of the Caputo derivative can be obtained 
by subtracting the singularity from the generator of the process obtained from an initial 
one by restricting the jumps to a chosen domain D. For instance, for a Feller process 
generated by the operator 

Af{x)=[ [f{x +z) - f{x) - {f'{x),z)x{z)]u{x,dz), (75) 


the process with jumps restricted to land on D (stopped-on-crossing-the-boundary 
process) for D the half-space fl53l) has the generator 


r*b—Xl 


Anfix) = 

k f dy2 f 




dy2 I dyiu{x-,y)[f{x+ y) - f{x) - {f'{x),y)x{y)] 

b — Xi 




dyiu{x-,y)[f {b,X 2 + 


' b—xi 


Vi 


2/2 - fix) - ifix),y)xiy)], (76) 


ioY X E D, which is a direct extension of fl5TD . However, unlike fl5T]l . A^fix) now diverges 
as X tends to a boundary point if df /dxiix) 7 ^ 0 there. In fact, as one sees directly, the 
main term of Ad fix) as ixi,X 2 ) ib,X 2 ) is 


dxi 


ixi,X2) 


dyi dy 2 ib - xi - yi)izib,X 2 ;yi,y 2 )- 


' b—xi 


/R'* 


Hence the analog of the Caputo derivative is 

df f 

ADi.fix) = ADfix)-—ib,X 2 ) dyi dya (& - a;i - |/i)z/( 6 , xs; 2 / 1 , 2 / 2 ). (77) 


dxi 


' b—xi 




4 Basic well-posedness results and examples 

4.1 The case /3 G (0,1), d = 1 

Let us start with the operator fl32|) assuming for simplicity that the kernels u have densities 
with respect to Lebesgue measure. Therefore let 

Afix)=xix)f'ix)+ [ ifix + y)-fix))uix,y)dy. (78) 

JR 
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Choosing for definiteness extension fl33|) (of (l3^ defined for x G (a, 6)) let 


A[a,b]*f{x) = 'y{x)f\x) + l^<b / (/[(x + y)Ab]- f{x))u{x, dy) 


“1“ [ {f[{x + y)y a]-f{x))u{x,dy), 

J—oo 


(79) 


which can be represented as the sum of the generators of a decreasing process, an increas¬ 
ing process and a drift: 


A 


[a,6]-*- 


Aa+i, + ^ 6 -* 


+ 7(0 


d 

dx 


with 


Aa+^f{x) = l 3 :>a / {f[{x + y) V tt] - f{x))u{x, dy), 

J —OO 


Ah_^f{x) = l^<b / {f[{x + y) Ab]- f{x))u{x, dy). 

Jo 

Let —oo < a < b < oo. Let us denote, as usual, by C[a,b] (resp. Coo(R)) the 
Banach space of continuous functions on [a,b] (resp. on R vanishing at inhnity), by 
Coo(—oo, a] (resp. Coo [a, cxo)) the Banach space of continuous functions on (—oo, a] (resp. 
[a, oo)) vanishing at inhnity, by C^(R), C^(—oo,a], C^[a, oo), C^[a,b] the subspaces 
of functions of the corresponding spaces Coo(R), Coo(—oo,a], Coo[a, oo), C[a,b], having 
derivatives up to order k from Coo(R), Coo(—oo,a], Coo[a, oo), C[a,6], respectively. 

Recall also that, for a domain D G R'^ with boundary dD and a Markov process Xx{t) 
in or just in D, a point x G dD is called regular if td^x) —)■ 0 in probability for 
X ^ Xq, X & D, where T£, is the exit time from D. We say that it is regular in expectation 
if EiTd{x) — 0 for X — )■ Xq, a; G D. 


Remark 5. 1. A point zero for D = (0, oo) is regular for a Brownian motion, but not 
regular in expectation, as F,td{x) = oo for all x & D. 2. We are using the notion of a 
regular point arising from the theory of parabolic PDFs (see e. g. or m), which is 
different from the corresponding notion used in the theory of Levy processes (see JE^). 


Theorem 4.1. Assume that'y{x) is continuously differentiable with a hounded derivative 
and that u{x, y) is a continuous function of two variables, which is continuously differen¬ 
tiable with respect to the first variable and has the following uniform bounds and tightness 
property 


sup 

X 


\y\^{x,y) dy < oo. 


sup 

X 



Au(x.y) 


dy < oo. 


(80) 


and 


limsup / \y\i'{x,y) dy = t). 

< 5^0 X J\y\<5 


(81) 


(i) Then the operator Aa+i, generates a Feller process on [a, oo) and a Feller semigroup on 
Coo[a, oo) with the invariant core C^[a, oo), the operator Ah_y, generates a Feller process 
on (—cx),6] and a Feller semigroup on Coo(—oo,6] with the invariant core C^(—oo,6]. 

(a) Moreover, the operator generates a Feller process on R and a Feller semi¬ 

group on Coo(R) and, ifjix) = 0, then A]^a,h]* generates also a Feller process on [a,b] and 
a Feller semigroup on C[a,h\ with the invariant core C^[a,h\. 
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(%%%) If 


or 


mm{\y\,e)u{a,y) dy > Cd’ 

O 

\m\\(y,e)v(b,y)dy > Ce’' 


(82) 

(83) 


for some C* > 0, r G (0,1), then the point a (resp. b) is regular in expectation for the first 
(resp. second) process in (i). 

(iv) If ^{ 0 ) < 0 (resp. 7 ( 6 ) > 0), then then the point a (resp. b) is a regular in 
expectation boundary point of the interval (a, b) for the process generated by 

(v) If A[a,b]* is the operator on the l.h.s. of fHB and > 0 , 7 jg > 0 , where (Ij^ is the 
unigue maximum of all (3j, then the points a and b are regular in expectation boundary 
points of (a, b) for the process generated by 


Proof, (i) Notice that the operator A]^a,b]* is a bounded operator C^[a^b] —)■ C[a,b] such 
that 

lim Aa+y,f{x) = 0, lim Ai,_^f{x) = 0 (84) 

x—>-a+ x^b- 

for / G C^[a, 6 ], as follows from 


lim (5 / iy{x,y)dy = 0. (85) 

< 5^0 J-R\[-5,S] 

Remark 6 . Eguation fl85|l is a conseguence of the first bound in flHOj) . In fact, since 

/ y^{x,y)dy = SF,,{6)+ F^{y)dy, 

Js Js 

where F^iy) = ^{x, z) dz, it follows that the both terms on the r.h.s. of this eguation 

are uniformly bounded. Hence, the l.h.s. of this eguation and the second term on the r.h.s. 
converge to yiy(x,u) dy implying 5Fx{5) —)■ 0 as 5 —)■ 0 . 

We shall follow the strategy of proof from Theorem 5.1.1 of [19]. Let us now work 
for dehniteness with A;,-* (other cases are dealt with analogously). Differentiating (and 
using straightforward cancelations) yields 


d 

dx 


Ab-y,f{x) 


[fix + y)- fix)] iy{x, y) dy - fix) / iy{x, y) dy 


' b—x 


r*b—X 


du 


du 

dx 


+ j ifix + y)- fix)) — ix, y) dy + (/( 6 ) - fix)) ^ 

Thus if / solves the equation 

/ = Ab-y^f, 

then g = f (if exists) solves the equation 

pb—x poo 

g= [gix+ y) - gix)]uix,y)dy - gix) / uix,y)dy 
Jo J b—x 


x,y) dy. 
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rb-x n rb /.OO 

+ dy g{x + z)dz—{x,y)+ g{z)dz —{x,y)dy. (86) 

t/ 0 V 0 k) X •) h—X 

Let us introduce the approximation Ab-i,h for our operator obtained by changing 
iz^x^y) to Vh^x^y) = l\y\^hi'{x,y). For any h the operator Ab_^h is bounded in C'(—cxo, 6 ] 
and hence generates a conservative Feller semigroup T/* there. Moreover, the operator on 
the r.h.s. of fl 86 l) becomes also bounded when u is replaced by so that this equation 
becomes well-posed in C(—oo, 6 ], so that is also a strongly continuous semigroup in 
C^(—cxo, b]. 

The key observation now is that T/* are bounded in C'^(—oo, b] uniformly in h, because 
the hrst two terms on the r.h.s. of (l 86 |l represent a conditionally positive operator with 
a negative coefficients at g, which therefore generates a positivity preserving contraction 
in C{—oo,b], and the last two terms are uniformly (in h) bounded operators. Hence 
{T^f)'{x) (by ' we denote the derivative with respect to x) are uniformly bounded for all 
h G (0,1] and t from any compact interval whenever / G C'^(—cxo, b]. Therefore, writing 

(ry - T^)/ = f‘TtUAi.-.h, - ds 

Jo 

for arbitrary hi > ^2 and estimating 


Jh 2 <\y\<hi 


<[ \\iTs"fy\W{x,y)\y\dy = o{l)\\f\\ci-oc,bh hi^O, 

Jo 

yields 

||(J,4. _ ft, ^ 0 . (87) 

Therefore the family T^f converges to a family Ttf, as h —)■ 0, which also forms a strongly 
continuous semigroup in C{—oo,b]. Writing 

Ttf-f Tif-T^f T^f-f 

t t t 

and noting that by fl87ll the first term is of order o(l)||/||ci(-oo, 6 ] as h —)■ 0 allows one to 
conclude that C^(—cxo, 6 ] belongs to the domain of the generator of the semigroup Tt in 
C(—oo, b]. 

Applying to Tt the procedure applied above to (differentiating the evolution equa¬ 
tion with respect to x) shows that Tt defines also a strongly continuous semigroup in 
C^(— 00 , 6 ], and hence C^(—oo, 6 ] is an invariant core for Tt. 

Remark 7. Notice that the semigroup Tt extends also to the strongly continuous semigroup 
on C'oo(R^ generated by with an invariant domain because the right condition 
of fl 8 T|) ensures a smooth gluing with the value A^-i^fx) = 0 for x > b for any f G C'^(R). 

(ii) This is quite similar and is omitted. 

(iii) To prove regularity of the boundary, we shall use the method of Lyapunov func¬ 
tions. Namely, to show that, say, 6 is regular for the last process in (i), it is sufficient to 
find a continuous function / in a neighborhood of [a, 6 ] such that / is differentiable for 
X G (a, 6 ), /( 6 ) = 0, and for x G (c, 6 ) with some c G (a, 6 ) one has f{x) > 0, A[a,b]*f{x) < 0 
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(see Proposition 6.3.2 of [E]). As such function one can take fu}{x) = {b — xY with some 
oj G (0,1). In fact, clearly fu^ip) = 0, f^x) > 0 for a: < 6 and for x approaching b from 
the left, Ai^a,b]*f{x) is of order 


-{b-x^-^ 



x)u{b,y) dy 


which tends to —oo as x ^ b for sufficiently small u under the assumptions of (iii). 

(iv) Using the same Lyapunov function fuj{x) = {b — xY one sees that the drift term 
always dominates the jump part part of the generator. 

(v) Proving regularity in a general case can be subtle. However, for the process from 
(v) the same Lyapunov function fcj(x) = (b — xY yields the required result. 

□ 


Theorem lO allows one to solve equations involving A\a,b]* or Aa±* by the standard 
techniques of stochastic analysis. Namely, let us denote by X^it) the Markov process 
generated by A of (1781) and by X*(f)[a, 6 ] the Markov process on [a, 6 ] generated by 
A^a,b]-k- 

Let us stress that the process X* (t) [a, b] (which is interrupted but not stopped at the 
boundary) is generated by A]^afi]* dehned on the domain C^(R) (or else on the domain 
C^[a, 6 ] if 7 ( 6 ) < 0 and 7 (a) > 0 ); the process X*(t)[a, 6 ; stop] on [a,b] obtained from 
X*{t)[a,b] by stopping at the boundary is generated by dehned on the domain, 

which is a subspace of C^[a,b] of functions / such that A[abp/(a) = A[abp/( 6 ) = 0 ; the 
process X* (t) [a, b; kill] killed at the boundary is generated by A\^a,b]-k dehned on the domain, 
which is a subspace of (7^[a, 6] of functions / such that A[ap]*/(a) = A[(j_ 6 ]*/( 6 ) = 0 and 
/(a) = fib) = 0. It is easily seen that if Jf^i'{a,y)dy = 00 and h'{b,y)dy = 00 , the 
operator Appp generates a strongly continuous semigroup (of the process X* (t) [a, b] killed 
on the boundary) on the subspace Co[a, &] of C[a,b] consisting of functions vanishing at 
a and b. 

Let T denote the hrst exit time for X*(f)[a, b] or X^it) from (a, b): 

T = inf{f > 0 : X*{t)[a,b] ^ (o, &)} = inf{f > 0 : X^t) ^ Y,b)}- 

Applying Dynkin’s martingale to a function / G C^[a,b] and Doob’s optional sampling 
theorem for the stopping time r (see e. g. [ 8 ] for the presentation of these two basic tools 
of stochastic analysis) yields 


f{x)=E f{X:{T)[a,b])- / {A[,,b]J){X:{s)[a,b])ds 

JO 

Hence if / is a solution to problem flT3|) . then 

/(x) =/(a)P(X*(r)[a,6] = a)+ /(6)P(X*(r)[a,6] = 6)-E f c/(X*(s)[a, 6]) ds. (88) 

Jo 

Moreover, since the trajectories of X^it) and X*(f)[a, 6] coincide till time r this can be 
expressed entirely in terms of the process Xx{t) as 


fix) = fiYPiXYr) <a) + /(6)P(X,(r) >b)-E giX^s)) ds. 


(89) 
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Introducing, as is common in the theory of Levy processes, the occupation-till-exit mea¬ 
sure H{x,dy) on (a, 6) by 


H{x,B) = e[ 1b 
J o 


{X^{s))ds, 


allows one to rewrite (j8^ as 


fix) = /(a)P(X.(r) < a) + /(6)P(X,,(r) > b) 



giy)Hix,dy). 


(90) 


(91) 


Thus we obtain the following result. 

Theorem 4.2. (i) Under the assumptions of Theorem \4.1\ problem 0431) can have at most 
one classical solution. //(x) < oo for all x E D, the probabilities P(Xa;(r) < a), 
P(X 3 ;(r) > b) are functions from C^[a, b] and the measure Hix, dy) is continuously weakly 
differentiable in x (so that the integral on the r.h.s. of fl^ belongs to C^[a, 6] for any 
g G C[a,b]), then formula (|^ supplies the unique classical solution to problem (H3l) . 

Remark 8. In the classical analysis of boundary-value problems for partial differential 
equations, problems like the one in diaD are usually understood to mean that the main 
equation there holds for all x excluding the boundary. The necessity of this attitude is 
easily seen here. Namely, if f belongs to the domain of the generator of a stopped or killed 
process, then the value of this generator on f at a boundary point should vanish. Thus 
only for g vanishing on the boundary the solution to 043 p can belong to the domain of the 
generator and satisfy the main equation up to the boundary. 


Furthermore, as is known from stochastic analysis, formula fl?Il) makes sense as a 
generalized solution under more general assumptions. Not going into much detail, let us 
only menton one particular situation. Namely, one says that a continuous function f{x) 
on [a, b] is a generalized solution to problem fl43P with = 0, if / belongs to the domain of 
the generator of the semigroup b] of the stopped process X*(t)[a, b; stop] on [a, b] 

(obtained by the closure of the operator A]^a,b]* dehned on the domain, which is a subspace 
of C^[a,b] of functions / such that A[a,b]*f{cL) = A[a,b]*f{b) = 0) and satishes A[a,b]*f = 0, 
or equivalently T/*°^[a, b]f = f. The following fact is a consequence of a general theory of 
boundary points (here the regularity of the boundary is crucial), see e. g. Theorem 6.2.3 
of [IS] for detail. 


Theorem 4.3. Under the assumptions of Theorem \4-l\ formula (ISip supplies a unique 
generalized solution to problem 


Furthermore, to solve problem fl44p one utilizes the process 


/(A'*(T)|a,6])e-"‘ + / e-"*(A - /l|„,,|.)/(A-(s)|a, 6]) <is, 

Jo 

which is known (see e. g. [S]) to be a martingale for any / G C~^[a, b] (under the conclusions 
of Theorem 14.11) . Again by the optional sampling theorem it follows that if / solves flT51) . 
then 


/W = E|/(.Y;(T)|a,(,])e-y 




(92) 
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Formula (l92|) again enures uniqueness of classical solution to and yields its integral 
representation in case the expectations involved in the r.h.s. of fl92|) are sufficiently regular 
functions of x. 

As an example, let us consider problem fH2]) on the interval [a,b] = [—1,1], The 
corresponding process Xx{t) of Theorem 14.21 is a symmetric Levy motion on R with index 
{3 G (0,1). For this process all ingredients of formula flTT]) are known, see e. g. [5]: 

P(A',(r) > 1) = 2 *-'’Iji _ (93) 

and H{x,dy) has the density 

where 

z = {l-x^){l-y^)/{x-yf. 

Thus fl?I|) yields a solution to problem in closed form. Some formulas for exit 
probabilities are also available for nonsymmetric Levy motions, see [3T] and [23], thus 
yielding explicit solutions to a slightly more general (compared to 042 p i problem 

aiD^+J{x) + a 2 D^_^f{x) = -g{x), /(a) = fa, /(6) = fb, (95) 


with arbitrary positive constants ai,a 2 - 

Similar results hold for the equations on a half-line involving the operators Aa±i,. 

4.2 The case /3 £ (0,1), (i > 1 

Let us now consider operator fl4Bp assuming again for simplicity that the kernel z/ has a 
density, h'{x,y), with respect to Lebesgue measure. 

The killed processes generated by (which are well studied for Levy processes, see e. 
g- 0) are generally easier for analysis than stopped processes. Therefore we concentrate 
on the analysis of operator (l50p . which is more involved. Let us consider only the case 
when D is the semi-space Db or the band D(^a,b), see (EH]), fl55l) . where is given by fl54ll 
and fl56ll . 

In what follows we have to use a rather ugly additional condition 

n{€,x) = e dyi dy 2 l\y\<iu{x;yi,y 2 )^ < Cu{e,x) (96) 

Je JRd Vl 

with a constant C, where 


u{e,x)= / dyi / dy 2 iy{x;yi,y 2 ). 
Je Jnd 


Its reasonability relies on the fact that it holds for stable-like processes in with 


p{x-,y) 


a{x) 

||^|cZ+l+/3 ’ 


P e (0,1), 


(97) 


20 






and hence for a variety of standard examples. 
In fact, since 




dy2'^\y\<l 


dr 


^arccos(r) 


sin'^ ^ dn, 


where dn is Lebesgne measnre on the nnit sphere in R'^ (or jnst a coefficient 2 in case 
d = 1) with the total area one has, assnming (ETj), that 


x) < a{x)\S'^ ^\e 


r 


dr 



d(p 


sin cj) 
cos^ 0 


and 


= a{x)\S'^ ^\e 



ajx) -/3 

l + fd 


uj{e, x) > a(x)|S'' 


d-li 


dr 


'0 


= a(x)|S"^ ^\e r ^ ^ dr dz = a{x) 


r‘arccos(r) 

sin (f) d(f) 
e-/3 _ 1 1 _ 


P 1-/3 


implying (l96D . 

Extending one-dimensional notations for fnnction spaces nsed above we shall denote 
by Coo [5] fhe Banach space of continnons functions on D vanishing at inhnity and by 
CUJD] its snbspace of fnnctions with Erst order partial derivatives belonging to Coo[D]. 


Theorem 4.4. Assume that n^x^y) = n{xi,X 2 ;yi,y 2 ) is a continuous function, which is 
continuously differentiable with respect to x and has uniform bounds flHOj) and tightness 
property flHTl) . where the integrals are over'Rf^^, and additionally the bound on the second 
derivative with respect to X 2 : 


snp 

X 



dxl 


^{xi,X2;yi,y2) 


dy < 00 . 


(98) 


Moreover, assume that fl96l) holds with a constant C and that n{x;yi,y 2 ) is an even 
function of y 2 - 

(i) Then, for a semi-space D = Db, the operator generates a Feller process Xf{x) 
on D and a Feller semigroup on Coo[.D] with invariant core CffD], 

(a) If additionally fl95]l holds also with the integral in yi taken over (— 00 , e), then also 
for a band D = 41(a,b) the operator A^^, generates a Feller process Xf[x) on D and a 
Feller semigroup on Ccc[D] with invariant core C^[Z)]. 

Proof. The proof follows the same lines as in Theorem 14.11 above. Let ns deal only 
with domain Db. Approximating o by Vh^Fy) = Iy 2 >h^{x]y) we get semigronps Th 
on Coo{D), which are nniformly bonnded and preserves twice continnons differentiability 
with respect to the second variable X 2 and bonnds to these derivatives. This holds becanse 
differentiation of the eqnation / = Ad*/ with respect to X 2 does not feel the bonndary 
so-to-say, that is we get 


d df _ r 

dt dx2 J^d 



dyiv{x]y)[ 


dx2 


{x + y) 


dx2 


(^)] 
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+ [ dy2 [ dyiiy{x;y)[-^ [b,X2 + 


/R^ 


' b—xi 


rb—xi 


+ 




dy2 


dX2 

du 


b — xi 
yi 


-y2 


dX2 




dyi^{x]y)[f{x + y) - f{x)] 
UX2 


/r^ 


dy2 


du 


' b—xi 


dyi-^{x]y)[f {b,x2 + 


b — xi 
yi 


1/2 -fix)], 


and similarly for the second derivative in X 2 (with h and withont it). The problem arises 
when differentiating the eqnation / = Ad*/ with respect to xi yielding the eqnation 


d 

df^ ~ 


b—xi 




dy2 


dyiu{x;y)[gix + y) - gix)]-g{x) dy 2 dyiu{x;y) 


/ r '* Jb-Xi yi ^^2 




u b-xi 

b, X2 H-2/2 

2/1 


' b—xi 


(99) 


(other terms cancel as in one-dimensional case) for g = df/dxi. Similar eqnation hods 
for Uh instead of u. Becanse of assnmed symmetry of u this rewrites as 


d 


f*b—xi 


dt 

with 


-gjgix) = / dy2 


dyiu{x;y)[g{x + y) - gix)] - gix)uib - xi, x) + (j){x; f) (100) 


fix;!) = 


'R'* 


a;(e,x) = / dy 2 dyiu{x;y), 

jRrf Je 


' b—xi 


yi 


dxl 


b, X2 + 6- — —y2 
2/1 


where 9 G (0,1) so that, for f jdx^ bonnded by a constant c, f is bonnded: 


ix; f)\ < cib - xi) / dy 2 / dyii 2 ix;y) 

jRd Jb-xi 


yi 

yl 


Eqnation fllOOp rewrites in the mild form as 

gt(x) = g„{x) + f‘ ds + f‘ •'>(‘-*>,^( 1 ; /) ds, 


( 101 ) 


where 


b—xi 


Agix) = / dy 2 
jR'i 


dyiuix;y)[gix + y) - gix)]. 


By fl96D . the last term in fllOip is nniformly bonnded, and hence eqnation fllOip and its 
versions with Uh instead of u have nniformly bonnded solntions for bonnded gQ. Hence we 
can now complete the proof as in Theorem 14.11 □ 

One can now get a direct mnlti-dimensional version of Theorem 14.21 for the bonndary 
valne problems 

Ad*! = 9^ UdD = 0 , ( 102 ) 

with g in D and f on dD given, which represent the simplest mnltidimensional analogs 
of linear eqnations with the Capnto derivatives. Alternatively, one can also analyze snch 
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problems via the reduction to killed processes (that is, to the analogs of RL derivatives), 
see Remark [U 

We shall not go into detail of general domains D here, but note that the problem 

AD^f = Xf + g, /|9^ = 0, (103) 

for A generating a stable Levy process in R'^ and D the ball H = {?/ G R'^ : \y\ < r} can 
be solved explicitly, using multidimensional extensions of formulas 0931) and 0941) given 
also in [5]. 


4.3 The case /3 G (1, 2) 

For the case (3 G (0,1) above we constructed the interrupted process on its own and then 
look at its stopping, which is quite natural. However, as we noted, the boundary-value 
problems (at least in one-dimensional case) for corresponding operators can be expressed 
in terms of the initial process stopped at the boundary. We shall follow this approach 
here, as the study of interrupted process becomes rather subtle. 

Let us reduce our attention to one-dimensional processes only generated by the oper¬ 
ators 

^fix) = J[f{x + y- f{x) - yf{x)]iy{x, y) dy 
with the density z/ satisfying 


sup 

X 


i\y\ A \y\‘^Mx,y) dy < co. 


'R 


(104) 


The question of whether such an operator generates a uniquely dehned process is non¬ 
trivial already in the case of this simple A, which can be looked at as the fully mixed-order 
fractional derivative. To go ahead, we shall use additional assumptions of regularity and 
monotonicity. The following statement is a particular case of Theorem 4.1 of mi: 


Proposition 4.1. Assume that v is twice continuously differentiable with respect to the 
first variable satisfying 

f d f 

sup / (li/l A l^ni—z/(a;,|/)| dy < oo, sup / (|y| A lyHl—z/(a;,y)| dy < cx), (105) 


and that the functions 


^{x,y), dy 


v{x,y) dy 


(106) 


are non-decreasing and non-increasing respectively for any a > 0. Then the operator A 
generates a Feller process Xfix) on R and a Feller semigroup with the space C'^(R) being 
an invariant core. The process Xfix) is stochastically monotone (but we will not use this 
latter fact). 


Next let —cxo <a<x<b<oo and let 

pb—x 


A[a,b]fix) = / [/(a V [(x + y) A 6]) - f{x) - y/'(x)]z/(x, y) dy 


(107) 
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be the operator representing the corresponding process X* (t) [a, b] interrupted on an at¬ 
tempt to cross the boundary of [a, b] (the processes on (— 00 , b] or [a, 00 ) with a one-sided 
boundary are considered analogously and will not be looked at) and 


Aia,b]*f{x) = A[a,b]f{x)-f{b) / [{b-x)-z]u{x,z)dz-f{a) 

J b—x 


[{a — x) —z]i'{x, z) dz 

(108) 


the corresponding analog of Caputo’s derivative (see d7T]) b 

As above, Proposition 14.11 allows us to apply the standard tools of stochastic calculus. 
Namely, let / G (^^[a,^] such that f'{a) = f'{b) = 0. Then we can continue it to all R by 
setting /(x) = f{b) for x > 6 and /(x) = /(a) for x < a and it will become a bounded 
continuously differentiable function on R. Denoting as above by r the exit time from 
(a, 6), that is r = inf{t : Xt{x) ^ (a, &)} and applying to / Dynkin’s martingale we get 
again (| 8 ^ . or else (jHI]), using the kernel H{x,dy) defined by (1^ and assuming / solves 
problem ([73]) with A = 0. 


Remark 9. Actually we can use Dynkin’s martingale only for twice continuously differen¬ 
tiable functions, and our (extended) f may have discontinuities of the second derivatives 
on the boundary, but this can be settled via approximation, as the final expression dUD 
does not involve the second derivative of f on the boundary. 


Therefore we get the following version of Theorem 14.21 for the present case jd G (1, 2): 

Theorem 4.5. (i) Under the assumptions of Proposition 4-l\ problem fl73|) with A = 0 
can have at most one classical solution. If the probabilities P(X 3 ;(r) < a), P(X 2 :(r) > b) 
are functions from C^[a, b] and the measure H{x, dy) is continuously weakly differentiable 
in X (so that the integral on the r.h.s. of belongs to C‘^[a,b] for any g G C[a,b]), 
formula fl9l|) supplies the unigue classical solution to fl73|) with A = 0. 

(a) Generally under the assumptions of Proposition \4. 1\ formula fl9l]) supplies a unigue 
generalized solution to fl73|) with A = 0. 


As was noted solutions to (Cl solve also flTijl under the additional assumption that 
the first derivative of the solution vanishes on the boundary, which yields a rather tamed 
(and expected) non-uniqueness for ([Till . 


5 Appendix 


For completeness we deduce here the expressions ([5|), (d and (d, IfTOj) from the original 
definitions. 

For (3 G (0,1) and x > a integration by parts yields 


iDm = 


r(i-/?) 


(x — t) ^f{t) dt = 


r d 

1 - 

1 

1 

1 _ 

la dt 

1 - 

1 

1 _ 


f{t) dt 


-nX-h 


(x — a) 

F(l-/3) 


so that 
A 


/(«) 


/(a) 


F(l-/3) 


rx 

(x — tf h 

L 

1 

1 

1 _ 


T{1-ld){x-a)h ' F(l-/?) 


+ 


f'{t) dt, 


(x — t) ^f'(t) dt. (109) 
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Another integration by parts using 


fit) = ^if{t)-f{x)), 

yields 

^/3 . _ fja) _ fif-fjx) _ r fit) - fix) 

^ T{1 - f3){x - a)t^ T{1-(3){x-a)t^ T{l-f3)J^ {x - t)^+t^ 

which equals to the r.h.s. of ([5]). On the other hand, 

D^+J{x) = la+^fix) = t)-t^fit) dt, 

which differs from (11091) by fia){x — a)~^/T{1 — (3) yielding (E]), ([7]). 

For (3 G (1, 2) and x > a one has 




( 110 ) 


which rewrites as 

= y{ 2-P) ^^ ~ f^~ffif - /'(«)) + -(3) / “ f~^ifif - fif) dt. 

Another integration by parts using 


fit) - fix) = J^ifit) - fix) -it- x)f'{x)) 

yields 

Da+Jix) = y{ 2-13) ^^ ~ af~ffix) - fia)) 


ix-a) ffia)-f{x)-{a-x)f\x)) + :^^^ 

which equals the r.h.s. of (ITn|) . 

On the other hand, again for /3 G (1, 2) and x > a, 


r fjt) - fix) -jt- x)f{x) 
Ja ix - 


t 2 —P 

'a+ 


fix) = 


r (2 - p) J, 

which rewrites by integration by parts as 


ix-t)^ ^fifdt, 


t2—0 

'a+ 


fix) = 


f{a) {x - a) 


- a)2-/3 


+ 


" {x -1)^-^ 


V{2-(3) 2-(3 r(2-/3)A 2-(3 

and by yet another integration by parts as 


-fit) dt, 


f(a) (x — a)'^ ^ f'(a)(x — a)^ ^ (x — t)^ ^f'if , 

r(2-/?) 2-/3 +r(2-/3)(2-/3)(3-/?)+/ r(2-/?) (2-/3) (3-/3) 
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Consequently, 


L>f+/(a:) = = 


-/ 3 ./ X , , r (x-ty 




r(2-/3) 


dt. 


Comparing this with flllOp yields ([9]) and flTT|) . 

Similarly, for /? G (1, 2) and x < a one has by dehnition flTTj) that 


Dtj{x) = 


r( 2 -/?) 


{a-xY ^{f{x) - /'(a)) 


{t — xY ^f”{t)dt 

E 

1-/9 


T{2-^Y ' " r(2-/3) 

By integration by parts this rewrites as 

1 


(t-x)-^if\t)-fix))dt. 


Dljix) = 


r(2-/3) 


{x - a)^ ^{f{x) - f{a)) 


r(i-/9) 


{a-x) ^(/(a)-/(x)-(a-x)/'(x))- 


r(-/9) 


fjt) - f{x) - (t - x)f{x) 
, {t-xY+^ 


dt, 


which equals the r.h.s. of fl2^ . Similarly, (1?I|) is obtained. 
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